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DESCRIPTION 
Artifact correction 

5 

The present invention relates to the field of image processing, for example in medical 
applications. In particular, the present invention relates to a method of artifact correction 
in a data set of an object of interest, to data processing devices and to respective 
computer programs. 

10 

In a CT scanner with a polychromatic source of radiation, such as a polychromatic x-ray 
source, a polychromatic x-ray beam passes through matter and low-energy photons are 
absorbed as the linear attenuation coefficient generally decreases with energy. As a 
result, the beam gradually becomes harder, i.e. its mean energy increases. The harder 
15 the beam, the less it is further attenuated. Therefore, the total attenuation is no longer a 
linear function of the thickness of the absorber. This non-linear behavior generally leads 
to a reduction in the attenuation coefficient. In bone for example, beam-hardening can 
cause reductions of up to 10 %. Neglecting this effect in the reconstruction process 
leads to well-known beam-hardening artifacts, such as cupping, streaks and flairs. 

20 

The reconstruction of CT images can be done with iterative algorithms. These 
algorithms, in particular iterative maximum-likelihood (ML) reconstruction algorithms, 
introduce artifacts along lines of high gradients of absorption values. This is well- 
known from the expectation maximation maximum-likelihood algorithm (EM-ML 
25 algorithm). Theses artifacts however are not desired, since they lead to a reduction of 
image quality. 

It is an object of the present invention to provide for an improved image quality. 



30 



According to an exemplary embodiment of the present invention as set forth in claim 1, 
the above object may be solved by a method of artifact correction in a data set of an 
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object of interest, wherein an image of the object of interest is reconstructed on the basis 
of the data set and wherein a statistical weighing is performed during reconstruction of 
the image. 

5 Advantageously, this may allow for a correction of artifacts in an image on the basis of 
a statistical weighing during reconstruction of the image. By weighting weighing certain 
data in the data set statistically the image quality may be improved significantly. 

According to another exemplary embodiment of the present invention as set forth in 
10 claim 2, the data set is a projection data set acquired by means of a source of 

electromagnetic radiation generating a beam and by means of a radiation detector 
detecting the beam. 

Advantageously, this may allow for an artifact correction of CT images acquired in a 
15 CT system, which may be used, for example, in the field of medical imaging, and 
therefore for an improved image quality. 

Another exemplary embodiment of the present invention is set forth in claim 3, wherein 
the source of electromagnetic radiation is a polychromatic x-ray source. Furthermore, 
20 according to an aspect of this exemplary embodiment of the present invention, the 

source moves along a helical path around the object of interest and the beam has one of 
a cone beam geometry and a fan beam geometry. 

The application of polychromatic x-rays may be advantageous, since polychromatic x- 
25 rays are easy to generate and provide for a good image resolution. Furthermore, since 
the geometry of the CT scanner system may be of different designs, such as, for 
example, cone beam or fan-beam geometry, a method of an exemplary embodiment-of 
the present invention may be applied to a plurality of different scanner systems and may 
not be limited to pure CT scanner systems, but also applicable to positron emission 
30 tomography (PET) or single photon emission computed tomography (SPECT). 
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According to another exemplary embodiment of the present invention as set forth in 
claim 4, the reconstruction of the image is performed on the basis of an iterative 
algorithm comprising a plurality of update steps until an end criterion has been fulfilled. 

5 Advantageously, by performing the reconstruction of the image iteratively, 

mathematical reconstruction methods may be used which comprise transcendental 
equations which can not be solved exactly. Furthermore, an iterative approach may 
provide for improved image quality, since the iterative loop may be repeated a plurality 
of times until the end criterion is met. 

10 

According to another exemplary embodiment of the present invention as set forth in 
claim 5, the iterative algorithm is a maximum-likelihood algorithm and the 
reconstructed image has the highest likelihood. Furthermore, the weighing is performed 
in each update step of the plurality of update steps. 

15 

Advantageously, a maximum-likelihood algorithm allows for determination of an 
equation which defines an image with the highest likelihood and a weighing in each 
update step of the plurality of update steps allows for an efficient compensation or 
correction of artifacts. 

20 

Another exemplary embodiment of the present invention is set forth in claim 6, wherein 
a number of detected photons during acquisition of the data set is determined and the 
weighing is based on a statistical error of the number of detected photons. 

25 Thus, the statistical error of the photon counts is estimated and taken into account 

during the image reconstruction, which may result in an improved image quality with 
— less artifacts. — 

According to another exemplary embodiment of the present invention as set forth in 
30 claim 7, an update of an attenuation parameter |Lij n+1 is calculated from the attenuation 
parameter |ij n by 
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y Zj M <l i ^>d i e^ 



wherein di is a number of photons emitted by the source of radiation, l i} is a basis 
function of an i~th projection, li is a vector of basis functions l s of the i-th projection and 
< L t , ju > is an inner product, defined by <l n fi >= ^ . 

This iterative statistical uncertainty weighted algorithm for maximum-likelihood 
transmission tomography may allow for an effective removal of artifacts by weighing 
the update with the intrinsic statistical error of the measured photon counts. 

According to another exemplary embodiment of the present invention as set forth in 
claim 8, the reconstruction of the image is based on a sub-set of at least two projections 
of all acquired projections of the projection data set. 

Advantageously, according to this exemplary embodiment of the present invention, 
performing the statistical weighing on sub-sets of projections may lead to a fast 
convergence of the image with less artifacts. 

According to another exemplary embodiment of the present invention as set forth in 
claim 9, a data processing device is provided, which comprises a memory for storing a 
data set and a data processor for performing artifact correction in the data set of the 
object of interest, wherein the data processor is adapted for performing the following 
operation: loading the data set and reconstructing an image of_the object of interest on 
the basis of the data set, wherein a weighing is performed during reconstruction of the 
image and wherein the weighing is based on statistical considerations. 
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Advantageously, this may allow for an improved image quality of images processed by 
the data processor which comprise artifacts caused by, for example, high attenuating 
objects, such as metal. 

According to another exemplary embodiment of the present invention as set forth in 
claim 10, the reconstruction of the image is performed on the basis of an iterative 
algorithm comprising a plurality of update steps until an end criterion has been fulfilled, 
wherein the iterative algorithm is a maximum-likelihood algorithm and the 
reconstructed image has the highest likelihood. Furthermore, the weighing is performed 
in each update step of the plurality of update steps. 

Advantageously, the iterative approach may allow for an improved and fast 
convergence of the image with less artifacts. 

According to another exemplary embodiment of the present invention as set forth in 
claim 11, a CT scanner system is provided, comprising a memory for storing a data set 
of an object of interest and a data processor for performing artifact correction in the data 
set of the object of interest, wherein the data processor is adapted for performing the 
following operation: loading the data set and reconstructing an image of the object of 
interest on the basis of the data set, wherein a statistical weighing is performed during 
reconstruction of the image. 

Advantageously, this may allow for an improved CT image quality of for example a 
part of the body of a patient, wherein the acquired data is processed by the data 
processor which comprise artifacts caused by, for example, high attenuating objects, 
such as a prosthesis or other implants. 

The present invention relates also to a computer program which may, for example, be 
executed on a processor, such as an image processor. Such a computer program may be 
part of, for example, a CT scanner system The computer program, according to an 
exemplary embodiment of the present invention, is set forth in claim 12. The computer 
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program may be preferably loaded into working memories of data processors. The data 
processors are thus equipped to carry out exemplary embodiments of the methods of the 
present invention. The computer program may be written in any suitable programming 
language, for example in C++ and may be stored on a computer readable medium, such 
as a CD-ROM. Also, these computer programs may be available from a network, such 
as the Worldwide Web, from which they may be downloaded into image processing 
units or processors, or any suitable computers. 

It may be seen as the gist of an exemplary embodiment of the present invention that 
artifacts in an image are corrected on the basis of a statistical weighing during the 
reconstruction of the image. This may lead to an improved image quality with less 
artifacts. 

These and other aspects of the present invention will become apparent from and 
elucidated with reference to the embodiments described hereinafter. 

Exemplary embodiment of the present invention will be described in the following, with 
reference to the following drawings: 

Fig. 1 shows a simplified schematic representation of an embodiment of a computed 
tomography (CT) scanner according to the present invention. 

Fig. 2 shows reconstructed image slices of a head phantom using the maximum- 
likelihood algorithm (left) and the statistical uncertainty weighted maximum-likelihood 
algorithm according to an exemplary embodiment of the present invention (right). 

Fig. 3 shows a flow-chart of an-exemplary embodiment of a method of artifact 
correction according to the present invention. 
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Fig. 4 shows an exemplary embodiment of an image processing device according to the 
present invention, for executing an exemplary embodiment of a method in accordance 
with the present invention. 

Fig. 1 shows an exemplary embodiment of a CT (computed tomography) scanner 
system according to the present invention. With reference to this exemplary 
embodiment, the present invention will be described for the application in medical 
imaging. However, it should be noted that the present invention is not limited to the 
application in the field of medical imaging, but may be used in other applications such 
as baggage inspection to detect hazardous materials, such as explosives, in items of 
baggage or other industrial applications such as material testing. 

The scanner depicted in Fig. 1 is a cone beam CT scanner. The CT scanner depicted in 
Fig. 1 comprises a gantry 1, which is rotatable around a rotational axis 2. The gantry is 
driven by means of a motor 3. Reference numeral 4 designates a source of radiation 
such as an x-ray source, which, according to an aspect of the present invention, emits a 
polychromatic radiation. 

Reference numeral 5 designates an aperture system which forms the radiation beam 
emitted from the radiation source to a cone-shaped radiation beam 6. The cone beam 6 
is directed such that it penetrates an object of interest 7 arranged in the center of the 
gantry 1, i.e. in an examination region of the CT scanner, and impinges onto the 
detector 8. As may be taken from Fig. 1, the detector 8 is arranged on the gantry 1 
opposite to the source of radiation 4, such that the surface of the detector 8 is covered by 
the cone beam 6. The detector 8 depicted in Fig. 1 comprises a plurality of detector 
elements. 

During a scan of the object of interest 7, the source of radiation 4, the aperture system 5 
and the detector 8 are rotated along the gantry 1 in the direction indicated by arrow 16. 
For rotation of the gantry 1 with the source of radiation 4, the aperture system 5 and the 
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detector 8, the motor 3 is connected to a motor control unit 17, which is connected to a 
calculation unit 18. 

In Fig. 1, the object of interest is disposed on a conveyor belt 19. During the scan of the 
object of interest 7, while the gantry 1 rotates around the item of baggage 7, the 
conveyor belt 19 displaces the object of interest 7 along a direction parallel to the 
rotational axis 2 of the gantry L By this, the object of interest 7 is scanned along a 
helical scan path. The conveyor belt 19 may also be stopped during the scans to thereby 
measure signal slices. Instead of providing a conveyor belt 19, for example in medical 
applications where the object of interest 7 is a patient, a moveable table is used. 
However, it should be noted that in all of the described cases it is also possible to 
perform a circular scan, where there is no displacement in a direction parallel to the 
rotational axis 2, but only the rotation of the gantry around the rotational axis 2. 

The detector 8 is connected to calculation unit 18. The calculation unit 18 receives the 
detection result, i.e. the read-outs from the detector elements of the detector 8 and 
determines a scanning result on the basis of these read-outs. The detector elements of 
the detector 8 may be adapted to measure the attenuation caused to the cone beam 6 by 
the object of interest 7. Furthermore, the calculation unit 18 communicates with the 
motor control unit 17 in order to coordinate the movement of the gantry 1 with motor 3 
and 20 are with the conveyor belt 19. 

The calculation unit 18 may be adapted for reconstructing an image from read-outs of 
the detector 8. The image generated by the calculation unit 18 may be output to a 
display (not shown in Fig. 1) via an interface 22. 

—The calculation unit which may be realized by a data processor may also be adapted to 
perform an artifact correction in the image based on read-outs from the detector 
elements of the detector 8. According to an aspect of the present invention, this 
correction may be performed by reconstructing an image of the object of interest on the 
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basis of the read-outs, wherein a statistical weighing is performed during reconstruction 
of the image. 

Furthermore, the calculation unit for the data processing device may be adapted for 



image on the basis of an iterative algorithm comprising a plurality of update steps until 
an end criterion has been fulfilled, wherein the iterative algorithm is a maximum- 
likelihood algorithm and wherein the reconstructed image has the highest likelihood. 
Furthermore, the weighing is performed in each update step of the plurality of update 



Furthermore, as may be taken from Fig. 1, for example, the calculation unit 18 may be 
connected to a loudspeaker 21, for example to automatically output an alarm. 

15 A standard approach for a reconstruction of CT images are iterative algorithms, for 

example a convex maximum likelihood algorithm. The Maximum Likelihood Method is 
able to determine an equation which defines an image with the highest likelihood L(/u) 
by setting dL({ij)/djUj = 0, where ju is the vector of attenuation parameters jUj . Since 
this equation is transcendental, it cannot be solved exactly, but an iterative approach can 

20 lead to the solution. One of these iterative approaches is the "Convex Algorithm" 

described in K. Lange and J. A. Fessler, "Globally convergent algorithms for maximum 
a posteriori transmission tomography", IEEE Trans. Med. Imaging 4, 1430-1450 
(1995), which is hereby incorporated by reference. The update for one iterative step is 
found to be 



5 performing artifact correction in the data set by performing the reconstruction of the 



10 steps. 



25 



J^l 9 <l i . t f>d l e 




(equation 1) 



where d i and Y t are the emitted and observed number of photon counts, respectively, 
Z f is the vector of basis functions l tj of the i th projection, and is the inner 



- 10- 



PHDE040070 EPP 



product 1-jjLij . This update ensures an increasing likelihood in the neighborhood of 
the optimal point, thus it should lead to an optimal image. 

The reconstruction result of such an iterative algorithm is shown as image slice 31 in 
Fig. 2. Image slice 31 is an image slice of reconstructed projection data of a forbild head 
phantom. However, as the reconstruction of the head phantom with noise shows, the 
algorithm according to equation 1 introduces artifacts along lines of high gradients of 
absorption values. This is well known from the expectation maximum likelihood (EM- 
ML) algorithm. 

The artifacts in the reconstructed image slice 31 occur as streaks. As may be seen from 
image slice 31 in Fig. 2, without any correction of the artifacts, i.e. without any artifact 
correction according to an exemplary embodiment of the present invention, a strong 
streak pattern can be observed. 

On the other hand, by reconstructing the image slice according to an exemplary 
embodiment of the method according to the present invention, in which the subtractions 
in the updates are weighted with the intrinsic statistical error cr Y of the measured 
photon counts Y i during reconstruction of the image slice, the artifacts are compensated 
for and the streaks disappear: 

* ~*' + *' Y,h<l,M->d t e^> • (e<JUaQOn2> 

where cr y . = is the statistical uncertainty of a Poisson distribution. Each of the 
updates is hereby the result of a respective cycle of the iterative algorithm. 
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It should be noted, that the method according to an exemplary embodiment of the 
present invention also works for a reconstruction with sub-sets of all measured 
projections. 

5 The image slices shown in Fig. 2 are the result of a convex algorithm (image slice 31), 
which is a special form of a ML algorithm, and a statistical uncertainty weighted 
improved convex algorithm (image slice 32), both after 500 iterations and starting from 
a uniform image. The level is 35 Houndsfield Units (HU), the window is 100 HU. 

10 The convex maximum likelihood algorithm used for reconstructing the image slices of 
Fig. 2 is described in K. Lange and J.A. Fessler, "Globally convergent algorithms for 
maximum a posteriori transmission tomography", IEEE Trans. Med. Imaging 4, 1430- 
1450 (1995), which is hereby incorporated by reference. 

15 Fig. 3 shows a flow-chart of an exemplary embodiment of a method of artifact 

correction according to the present invention. After the start in step SI, an acquisition of 
the projection data set of the object of interest is performed in step S2. This may be 
performed by detecting a polychromatic x-ray beam emitted from a radiation source, 
which moves along a helical path around the object of interest. The x-ray beam passes 

20 the object of interest and is attenuated by object points inside the object of interest. 
During the movement of the radiation source a plurality of projections are acquired, 
each projection belonging to a respective mean position of the radiation source. After 
acquisition of the projection data set of the object of interest the number of detected 
photons Yi of an i-fh projection is determined in step S3, e.g. with the help of a data 

25 processor. Then, in step S4, the statistical error <j y of the number of detected photons Yi 
is calculated by the data processor and the vector \i of attenuation parameters jij is set to 
an initial value, e.g. /u j = 1 , in step S5. Each attenuation parameter Mj belongs to a 
respective interval along the projection line of the i-th projection. 
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In step S6 then, a counter is set to n = 0 after which the iterative loop for the image 
reconstruction starts. In step S7 the n-th update (which is in case of the beginjoing of the 
iterative loop the first update) of the attenuation parameter jlij is calculated by 

Z^ru v , / 2 

2Li l ij<li>M n >d i e^ > 

Then, the method continues to step S8, where it is determined whether an end criterion 
is met or not. The end criterion may, for example, be that the iteration has been 
performed for a certain predetermined number of times or that a difference between 
consecutive updates does not exceed a certain predetermined threshold value. In case 
that it is determined in step S8 that the end criterion is not met, the counter is increased 
by one (step S9) and the method continues with step S7, where (n + l)-th update of the 
attenuation parameter jlij is calculated. Step S7 - S9 may be iteratively repeated until the 
end criterion is met. In case it is determined in step S8, that the end criterion is met, the 
method continues to step S10, where it ends. 

It should be noted that, as may be apparent to the skilled person, the above described 
technique may be applied to all known reconstruction techniques in which a difference 
is summarized such as in the convex algorithm described above. 

Fig. 4 depicts an exemplary embodiment of a data processing device according to the 
present invention for executing an exemplary embodiment of a method in accordance 
with the present invention. The data processing device depicted in Fig. 4 comprises a 
central processing unit (CPU) or image processor 151 connected to a memory 152 for 
storing an image depicting an object of interest, such as a patient. The data processor 
151 may be connected to a plurality of input/output network or diagnosis devices, such 
as an MR device or a CT device. The data processor may furthermore be connected to a 
display device 154, for example a computer monitor, for displaying information or an 
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image computed or adapted in the data processor 151. An operator or user may interact 
with the data processor 151 via a keyboard 155 and/or other output devices, which are 
not depicted in Fig. 4. 

5 Furthermore, via the bus system 153, it is also possible to connect the image processing 
and control processor 151 to, for example a motion monitor, which monitors a motion 
of the object of interest. In case, for example, a lung of a patient is imaged, the motion 
sensor may be an exhalation sensor. In case the heart is imaged, the motion sensor may 
be an electrocardiogram (ECG). 

10 
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CLAIMS 

1 . A method of artifact correction in a data set of an object of interest, the method 
comprising the step of: 

reconstructing an image of the object of interest on the basis of the data set; 
wherein a statistical weighing is performed during reconstruction of the image. 

2. The method according to claim 1, 

wherein the data set is a projection data set acquired by means of a source of 
electromagnetic radiation generating a beam and by means of a radiation detector 
detecting the beam. 

3. The method according to claim 2, 

wherein the source of electromagnetic radiation is a polychromatic x-ray source; 
wherein the source moves along a helical path around the object of interest; and 
wherein the beam has one of a cone beam geometry and a fan beam geometry. 

4. The method according to claim 1, 

wherein the reconstruction of the image is performed on the basis of an iterative 
algorithm comprising a plurality of update steps until an end criterion has been fulfilled. 

5. The method according to claim 4, 

wherein the iterative algorithm is a maximum likelihood algorithm; 

wherein the reconstructed image has the highest likelihood; and 

wherein the weighing is performed in each update step of the plurality of update steps. 

6. The method according- to claim 2, further comprising the step of: — 
determining a number of detected photons during acquisition of the data set; 
wherein the weighing is based on a statistical error of the number of detected photons. 
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7. The method according to claim 5, further comprising the step of: 
determining a number of detected photons Yi during acquisition of the data set; 
wherein the weighing is based on a statistical error <j y . of the number of detected 
photons Y } ; 

wherein an update of an attenuation parameter uj n+1 is calculated from the attenuation 
parameter by 

Zu i if i / 2 

m" +1 = Mj + M- Z^lLZlL 

wherein di is a number of photons emitted by the source of radiation; 
wherein ly is a basis function of an i-th projection; 
wherein li is a vector of basis functions ly of the i-th projection; and 
wherein <l i9 jn >= ^ is an inner product. 

8. The method according to claim 2, 

wherein the reconstruction of the image is based on a sub-set of at least two projections 
of all acquired projections of the projection data set. 

9. A data processing device, comprising: 

a memory for storing a data set of an object of interest; 

a data processor for performing artifact correction in the data set of the object of 
interest, wherein the data processor is adapted for performing the following operation: 
loading the data set; 

reconstructing an image of the object of interest on the basis of the data set; 
wherein a statistical weighing is performed during reconstruction of the image. 
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10. The data processing device according to claim 9, 

wherein the reconstruction of the image is performed on the basis of an iterative 
algorithm comprising a plurality of update steps until an end criterion has been fulfilled; 
wherein the iterative algorithm is a maximum likelihood algorithm; 
5 wherein the reconstructed image has the highest likelihood; and 

wherein the weighing is performed in each update step of the plurality of update steps. 

1 1. A CT scanner system, comprising: 

a memory for storing a data set of an object of interest; 
10 a data processor for performing artifact correction in the data set of the object of 

interest, wherein the data processor is adapted for performing the following operation: 
loading the data set; 

reconstructing an image of the object of interest on the basis of the data set; 
wherein a statistical weighing is performed during reconstruction of the image. 

15 

12. A computer program for performing artifact correction in a data set of an object of 
interest, wherein the computer program causes a processor to perform the following 
operation when the computer program is executed on the processor: 

loading the data set; 

20 reconstructing an image of the object of interest on the basis of the data set; 

wherein a statistical weighing is performed during reconstruction of the image. 



25 
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ABSTRACT 
Artifact correction 

5 The reconstruction of images of an object of interest may introduce artifacts along lines 
of high gradients of absorption values. According to an exemplary embodiment of the 
present invention, these artifacts may efficiently be removed by a statistical weighing 
during reconstruction of the image. Advantageously, according to an aspect of the 
present invention, the reconstruction of the image may be performed iteratively, 
10 wherein the updates are weighted with the intrinsic statistical error of the measured 
photon counts. This may lead to an efficient removal of artifacts. 
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